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Математическое моделирование процессов 
кристаллизации металла 


Исследуется пространственная задача Стефана с учетом конвективного движения и примесей в жидкой 
фазе. Получено уравнение свободной границы. 


1. Работа посвящена изучению процессов кристаллизации двухкомпонентных сред 
в случае, когда распространение тепла связано не только с теплопроводностью, но и 
с конвективным переносом, присутствующим в жидкой фазе вещества. Рассматривае- 
мая задача включает в себя как двухфазную задачу Стефана, так и начально-краевую 
задачу для системы Навье-Стокса, описывающую движение вязкой несжимаемой жид- 
кости в нецилиндрической области. При изучении задачи учитывается скачек плот- 
ности вещества на границе раздела фаз. 

Для описания поля скоростей в зоне поступления перегретого металла исполь- 
зуется математическая модель затопленной струи вязкой жидкости, основанная на 
известном в теоретической гидродинамике точном решении нелинейной системы диф- 
ференциальных уравнений Навье-Стокса, а поскольку в затопленной струе переме- 
шивание можно охарактеризовать так называемой свободной турбулентностью, ха- 
рактеризующейся одним числовым параметром — коэффициентом «кажущейся», или 
турбулентной, вязкости Г’, , полученные формулы можно интерпретировать как в ла- 


минарном, так и в турбулентном приближениях. 

Упомянутое решение системы Навье-Стокса имеет полярную особенность в на- 
чале координат и обладает осевой симметрией, а его функциональная структура на 
столько проста, что для соответствующего уравнения конвективной теплопередачи 
удается получить решение с двумя произвольными вещественными параметрами, оп- 
ределение которых приводит к полному описанию температурного поля в пределах 
исследуемой зоны. В сферической системе координат (хи, 6, ф, ) радиальная и азиму- 


тальная компоненты скорости выражаются по формулам 


и. а - 
гзшд 


у 
|: > 
зто $1 а+1—-с0$0 


и = 


в которых у — кинематическая вязкость среды, а параметр а однозначно выражается че- 
рез полный поток импульса Р рассматриваемой струйки, причем 


16лу? 552“ В? 
ВИА | если РОТ, ей. 
ты Р 
где р -— плотность жидкости. Полный поток импульса струи Р можно приближенно 


‚ если РОТ, 


посчитать по формуле 
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где С — количество поступающего расплавленного металла за одну секунду, а Г, — аб- 


солютное значение скорости движущегося вниз металла на оси кристаллизатора в точке 
входа струйки в металлическую ванночку. Температурное поле описывается формулой 


В 2а 
Т(".6)= | ы аа. 
г \а+1|-с039 


в которой о =урс/ А — число Прандтля, характеризующее свойство жидкой среды, а 


Ви С - вещественные параметры интегрирования, однозначно определяющиеся че- 
рез тепловую интенсивность О = ССТ, , нагретой до температуры Т, струйки и число 
Нуссельта @ =аК/А. 

Вычисления показывают, что 


2а 
о ВИ 
Ю\а-+1 ав 


где Г — температура внешней относительна кристаллизатора среды, а постоянная В 
в силу асимптотических соотношений определяется по формулам: 


о Ба О. если РОТ, 


4лрсу За(а+2) 


__ = В(1-а”), если РО1. 
8лосу 


Каждая линия тока рассматриваемого поля скоростей имеет единственную бли- 
жайшую к оси кристаллизатора точку и все такие точки лежат на одной и той же пря- 
мой, имеющей некоторый угол наклона 6, к оси симметрии, а вращение той прямой 
вокруг оси описывает конус — боковую границу рассматриваемой зоны. Другим харак- 
терным размером зоны является глубина 7, жидкой ванночки вдоль оси кристаллиза- 


тора. Как показывают непосредственные вычисления, 0, и 7, определяются из формул 


где Г” — температура кристаллизации. 

Для иллюстрации приведенных формул предположим, что высота шлаковой ван- 
ночки равна 10 см, а плотность шлака в 2,5 раза меньше плотности металла, и поскольку 
металл поступает в виде серии капель приблизительно 7 грамм (следовательно, и мас- 
са капли равна т, =7 г), то каждая из них, если учитывать только силы тяжести и ар- 
химедовы силы, в момент соприкосновения с зеркалом жидкого металла имеет ско- 
рость около 110 см/с и несет импульс Р, = 770 г см/с. Производительность печи в одну 


тонну веса в час соответствует С = 280 грамм веса в секунду, следовательно, в одну 

секунду выпадает около 40 капель, суммарный поток импульса которых равен прибли- 
4 2 

зительно Р = 3х10` гсм/с°. И для того случая, когда все капли сливаются в сплошную 


струйку вдоль оси кристаллизатора, несущую поток импульса Р = 3х10* г см/с?, и 
для того случая, когда в некоторую точку за каждую секунду опускается одна капля 


1 
с импульсом Р, = 770 г см/с, формула с0$6, = ЕЕ: дает нулевой угол 0, , следователь- 
а+ 


но, зона интенсивного перемешивания на значительной глубине жидкой металличес- 
кой ванночки имеет вид цилиндра, радиус которого, как показывают наблюдения за 
зоной интенсивного каплепадения, равен приблизительно одной трети радиуса электрода. 
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Если считать, что металл поступает перегретым на 150 °С, так что его температу- 
ра приблизительно равна Т., =1650° С, то по указанной выше формуле получаем теп- 


ловую интенсивность О =3,9х10* кал/с. Пользуясь формулами, получаем и. =7А, 
если при этом принять во внимание, что для расплавленных металлов молекулярная 
вязкость приблизительно равна у = 5х10`* см/с, а число Прандтля имеет близкое к 0,08 
значение. Применительно к ЭШП полученная величина имеет условный характер, 
прежде всего потому, что металл в этом случае поступает не в виде сплошной струйки, 
а в виде серии капель. Применяя предыдущие формулы и расчеты к случаю, когда в 
одну и ту же точку ежесекундно падает одна капля с указанными параметрами, полу- 
чим 7; = 0,17^А = 8, 5 см, то есть около одной пятой радиуса. Истинная глубина при ЭШП, 
несомненно, больше, поскольку мы не учли теплопередачи к жидкому металлу со сто- 
роны высокотемпературной шлаковой ванночки. 

Существенную поправку к указанным величинам вносит эффект турбулентности, 
который в рассматриваемом случае «свободной турбулентности» описывается коэф- 


фициентом «кажущейся» вязкости И, . Согласно известной формуле И, =11.10°/К, 


Р 
где К =-— - поток кинематического импульса. Для случая струйного течения получа- 
р 


ем И, = 0,215, т.е. Г, = 140 и, тогда, как в случае серии капель, весьма правдоподоб- 
ным является соотношение Г, = 10у=5х10`? см * /с. 


Переходя к турбулентному переносу тепла, примем в соображение, что соответст- 
вующий коэффициент вдвое превосходит коэффициент турбулентного переноса импульса. 
Уравнение энергии в турбулентном потоке в этом случае характеризуется «кажущимся 
коэффициентом теплопроводности» К = 27, + К ‚где К -— «коэффициент молекулярной 

и 
(К-+2И.) 


В разобранных выше случаях имеем о’, = 0,5 для струйки и о’, = 0,31 для серии капель, 


температуропроводности», а соответствующее число Прандтля равно о; = 


что в первом случае дает г, =3,5 см=7х10? А, аво втором 7; = 12ем=2,4х10? В. 


В условиях ЭШП эти величины увеличиваются из-за потока тепла со стороны ванноч- 
ки жидкого шлака. 


2. Пусть О, — заданная область в №*, граница которой состоит из двух замк- 
нутых связанных гладких поверхностей Г,’ и Г, , не имеющих самопересечений. 
Пусть, далее Г. — гладкая замкнутая поверхность, лежащая внутри <>», , такая, что Г. 
лежит внутри ограниченной области, границей которой является Г, . Поверхность Гу 
разбивает ©, на две подобласти О.’ и ©, , которые в начальный момент { = 0 заня- 


ты жидкой и твердой фазами соответственно. Будем обозначать через ©” область 


занятую жидкой (твердой) фазой в момент времени 1. Заметим, что в процессе кри- 
сталлизации проходит изменение границы Г.’ (это связано с тем, что жидкая и 


твердая фазы имеют разные плотности), а граница Г, остается неизменной. Задача 


состоит в определении областей О” и О’ (т.е. границ Г и Г, ), занимаемых твер- 
дой и жидкой фазами соответственно в момент времени Ее [0,7], вектора скорости 
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Т(х, 1 =(И(%,0),И, (х,1), Г, (х,1)) ‚ давления р(х,г) , концентрации примеси с(х,{) , распре- 


делений температур жидкой и’ (х,ё) и твердой и (х,г) фаз по следующим условиям: 


и 


+(ГУ)и* (ха? Уи (х,6 = 0, (х,1) Е 2%; —а?У‘и (х,6)=0, (х,)ЕРх, 


ди (х,г) 
ГА 
от Руб В+ Уре, Г) = №? Г(х, 0 Ги ‚с), УГ(х, 1) =0, (х,6) ЕД», 


Р(х,0) = Со: Т(Р, руп = -абьд, (х,) ЕГИЗИ, =-а- РТ. = 0, (0 еГ, 
р 


и*(х, 6) = В*(х, г), (х,6)ЕГ/ УГ, ; и*(х,0) = А*(х); (1) 
и’ =и и Е НИ Й’, (х,1)ЕГ,, 

ди ди 
дс(х,1) 


Я (РУАН Й- И с(х, 0, (х, 0) ЕД, *;с(х,0) = 2, (х), 


2 = абы дс) еГиз-а = Вет, „бьд ет, 
11 


Здесь О;^ = {(х,1):хЕ О, ЕЕ (0,Т)}, х=(х,х,,х,), О*, — области соответственно жид- 

кой и твердой фаз, дО* =Г‚чГ*, 60 =Гу Г, п - нормаль к Г, ‚ направленная в 

сторону О*; Т’, и, &, д, р’, р, а, В, у, к,к, — положительные постоянные, 
9" :0 ОУ. 

= ее —,—) Т (И, Р)- тензор напряжений с элементами Т. +—^), 


д, ей 
ОХ Ох. ’ дх, . дх, дх 


У 7 


И’ — скорость движения фронта кристаллизации в направлении нормали п, Ги! - 
нормальная и тангенциальная составляющие Г. Если Фо =и(хб+ес(х-Т’=0- 
У(и* + Ес) 

[Ум +=с) |. 

Укажем, что условие Стефана можно представить также в виде: к *|Уи [р - 
—к. | Ум [ +2=к *(Уи „Ус)-2=к.^(Уи*,Ус) = ур" (ки, +к.и’,), (х,ЕГ,. 

2+а 
Предполагается, что АС) ЕН““(©,”), С(х) с И?" (© *), В*(ьдЕН“" 2 (ГА ОГух 


2+а 


х[0,7]), 1 (м, с) Е СКВ”), в(хбеЕН”" 2 (Г*х[0,Т]), 8) ЕН“ (бу ). При этом 


уравнение поверхности Г, , тогда И’, =-—Ф, /|УФ |, п= 


1+а 


&(х,1) и &, (х,Г) должны быть функциями класса Н Ви (®*х[0,Т]). Предполагается 


также выполненными условия согласования до первого порядка включительно, кото- 
рые следуют из предположения существования гладкого решения и формулируются 
аналогично [1, с. 268, с. 363]. 


Отметим, что при малых значениях ([, задача (1) разрешима в классе В функ- 
2+а 2 2+а 


а 


2+а 


2+а, 
ницы Г’, и Г, описываются функциями, принадлежащими классам Н а 
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Решение задачи (1) моделирует процесс кристаллизации вещества с учетом пе- 
реноса примеси в жидкой фазе. При этом последнее условие в (1) следует из закона 


Нернста, а /(и’,с) описывает влияние неравномерного распределения температуры 
и концентрации примеси на движение жидкости. 


3. Известно, что свободные границы Г, и Г’ можно представить в виде Г, = 
= {х = х(@)+п(@)р(®,0}, Г’ ={х=х(9')+1(@',0)п(’) ‚ где @=(@,6,), @ =(@', 


©’, ), х(@)ЕГу, х(@)ЕГ,, р(@,1) и 1(@’,Г) — некоторые функции соответственно 
аа: 2+а, 2+4: ре х 
классов НН * (Г. х[р0,ТриН *? (Г, х 0,Т]), р(@,0) =0 и п(@ ,0)=0 [2]. 
Предложен метод решения задачи (1), состоящий в разложении решения в ряд 
по степеням малых чисел &: 


и* (х,1; =) Нид + У и (0) ‚РО, Е)= в) +У =‘ ри (к, , 


К= К=1 


Ив) =, +У ЕТ, (хи = 12,3; рав) =У =“ р, (0), (2) 
Е=1 


К=1 


с(х, 1) = с (х)+ у соо 


Для нулевого приближения и“, (%),И, (>) = (И, (<), Г (х), Г (х)), Го и с, (х) из ус- 
ловий (1) и разложения (2) вытекает следующая задача: 
(ИУ), (о) + Урь (о) = И (+ У (шь,сь), хе Оу, 


УР. (х)=0, хЕО,*, ТИ, ри =-9 (п, хеГу, И =-а- РМ, И =0, 


п р 
хеГь; (ИУ)и* —а,?Уи* =0, хеО, и (х) = В*(х), хеГ Го, (3) 
. ди . 
и, (Х=щ, (>) =Т, хЕГ,; К. — К. В 0, хЕГ., У*и, =0, хеО,, 
ди ди 
(ИУ) — с, =0, хеО,", с, (х) = в, (х), хеГ,; аа = ЕЁ. 
и 


Здесь ради простоты предполагается, что функции В* и 4 зависят только от пе- 
ременной х. 

Лемма 1. Пусть функции и,’ (х)= А*(х), И, (х)=С(%), с, (х)= 8, (х) являются 
решением задачи (3) соответственно в области 0’ и О’’. Тогда эти функции можно 


взять в качестве нулевого приближения задачи (1). 
4. Далее, пусть 


ог =О, [0,7], Г =Го х[0,Т], Ги =Го [0,7], Ги =Г.х[0,Т]. 


+ 
Рассмотрим первое приближение (И,и, ‚р, р,,с,) задачи (1) для малых чисел 5. 
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Имеем: 


ый +(ГУ), +(ГУ)И. + Ур, =? [А р, ‘(м сои“ р, (ис) (ВОт ь; 


УХ = 0, (х,6) Е О Е 


Т(,+Г, р)п=0,хЕГ,*,Й (х,0) =0, ее 
т р =0,хЕГ,; 
и. 4 (ГУ)и + (ИУ)и“ —а `У‘и! =0,(х,) ЕО)"; = 0, (х,1) ЕО, , 
и*(х,0) = 0;м (х,Е) = 0, (хЕГ Г” г, (5) 
и =и, К о 2%. нЕ р = + Д(х, В) = ур" =. (х,1) Г; 
е +(У)с, +(ИУ)с, - №? =0,(х,) ЕО", 
с (х,0) =0,с, (х,1) = 0,(х,1#) Е Г*х, 
наи = 2 (%,0, (хдеГу, КО)= аи . — 
и пе, +с, (х,6)=0, (Гу. 


Зададим теперь Г = Й,(х,В. Затем решим задачу (5), (6) и найдем и*,с, р. Пос- 
ле чего заменим и*,с, р — решением задачи (5), (6) и решим задачу (4), являющуюся 


начально-краевой задачей для системы Навье-Стокса. Затем, используя новое 
значение У (х,!), снова решаем задачу (5) и (6) и т.д. Таким образом, получим процесс 


последовательных приближений. Доказательство сходимости этого процесса анало- 
2+а 


гично приведенному в работе [3]. При этом при заданном р (4,1) ЕН в (Гог) най- 


2+а, +4, 
дем функции и*,(х,6р)ЕН °? (0*,), с (х,6р)ЕН ? (0*,), как единственное 
решение задачи (5), (6) [1], причем р’ (©, о находим как неподвижную точку сжимаю- 


2+а 2+а 


ди* х 5С 
щегося оператора М, : М,р, = а 2, и + Л(х)аьЬ х(@) Е Гуу. 
дп 2" т 
Имеют место следующие утверждения. 
Лемма 2. Пусть выполнено условие | УД" (х) |= ние наГ,. Тогда оператор М. , 
т 


2+а 2+а 


а.2+2 ‚2+а 
действующий из Н '° (Гу) вН * (Гу), имеет там неподвижную точку. 


Лемма 3. В качестве первого приближения задачи (1) можно взять решение за- 
дачи (4) — (6): и’ (хуйуе, 5В,Й 5 ),р, Сы, ю, (х/). 
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ТМ 


Теорема. Пусть Но 
т 


20 на Е. . Тогда при малых числах & и достаточно ма- 


лых значениях { справедливы формулы: 


т. = ш' (х(6),1) 
Г: х=х(@) ее ЕТит, 
ны * — < (х(@),0) + 8 (х(@))- 8(х(@),1) + 
Г”: х=х(@ )-=п ЕС +0(=),(х,1) Е Г, 
дп 


2+а 


2-+а, 
где и“ (х,1), с, (х,0, р, (@,0), 1, (9,1) — функции класса Н `? , являющиеся ре- 
шением задачи (4) — (6). 
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О.С. Мненко 

Математичне моделювання процесв кристалзацй металу 

Дослиджуеться просторова задача Стефана з урахуванням конвективних рух!в 1 домппок у рдинай 
фаз1. Доведено равняння в1льно! границу. 


А.5. Мтепко 
Маветайса|! Моде!по о! пе Ргосез$е$ СгузбаШтайоп о? Меёа1 


Тве @гее Аппепз1опа| сопуесноп З{еЁап рго ет ш Па рВазе 15 шуезиожщеа. Когли]а оР Йее Боипдагу 
едиайоп 1$ оМатеч. 
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